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LONG-TERM  GOALS  &  OBJECTIVES 

Many  owl  species  rely  on  specialized  plumage  to  achieve  functionally  silent  flight,  which  requires  the 
suppression  of  aerodynamic  noise  for  frequencies  2+kHz  to  not  affect  their  acoustical  targeting  system 
and  not  alert  prey  of  their  approach.  The  specialized  plumage  features  believed  to  be  responsible  for  this 
noise  suppression  are  a  leading-edge  comb  of  evenly  spaced  feathers,  the  velvety  down  material 
distributed  over  the  upper  wing  surface,  and  a  tattered  and  compliant  trailing  edge  (see  figure  1).  Recent 
theoretical  work  by  Jaworski  and  Peake  [1]  suggests  that  trailing-edge  noise  can  be  mitigated  by  the 
porous  and  elastic  properties  of  the  edge  to  the  point  where  turbulence  scattering  by  the  edge  scales 
similarly  to  roughness  noise.  If  validated  experimentally,  this  finding  has  broad  implications  for  the 
design  of  quieter  airframes  and  rotors,  and  the  diminished  edge  noise  implies  that  roughness  noise  has  a 
more  substantial  contribution  to  the  acoustic  signature  of  specially-designed  blades.  Thus,  the  focus  of  the 
present  research  is  to  determine  more  precisely  how  these  blades  must  be  designed  to  achieve  the  desired 
turbulence  noise  reductions,  and  to  model  and  assess  how  the  novel  surface  roughness  of  owls  may 
further  contribute  to  tunable  acoustic  stealth. 

In  contrast  to  unique  owl  features  at  the  leading-  and  trailing-edges,  the  velvety  down  is  the  least  studied 
and  most  poorly  understood  component,  which  according  to  Lillcy  [2]  is  the  feature  that  may  offer  the 
means  to  reduce  the  generation  of  turbulence  noise  at  its  source.  The  length  scale  of  the  owl  down 
(~lmm)  is  inconsistent  with  the  range  of  frequencies  for  a  conventional  acoustical  absorber,  thus  the 
noise  reduction  mechanism  is  likely  due  to  a  novel  coupled  interaction  between  the  flexible  downy  fibers 
with  the  fluid  flow.  In  addition  to  determining  theoretically  and  validating  experimentally  the  parameter 
space  where  trailing-edge  noise  scales  the  same  as  or  weaker  than  roughness  noise,  the  present  modeling 
effort  at  the  University  of  Cambridge  seeks  to  identify  the  salient  parameters  of  roughness  noise  due  to 
flexible  surface  elements.  This  effort  seeks  to  establish  an  understanding  of  how  the  noise  may  be 
mitigated  differently  when  the  down  is  modeled  either  as  a  collection  of  individual  flexible  fibers  or  as  a 
poroelastic  continuum. 

This  project  is  part  of  a  collaborative  effort  between  Lehigh  University,  University  of  Cambridge  (UK), 
Virginia  Tech,  and  Florida  Atlantic  University.  Incipient  funding  under  the  NICOP  program  (N62909-12- 
1-7116)  to  the  University  of  Cambridge  supported  the  theoretical  analysis  on  the  noise  generation  by 
flexible  roughness  by  Drs.  Nigel  Peake  and  Justin  Jaworski  (now  at  Lehigh  University).  Researchers  at 
Virginia  Tech  (Dr.  William  Devenport  and  graduate  student  Ian  Clark)  continue  to  carry  out  experimental 
work  to  test  the  developed  trailing-edge  and  hydrodynamic  flexible  roughness  theories  and  provide 
support  for  further  theoretical  developments.  Dr.  Stewart  Glegg  at  the  Florida  Atlantic  University  advises 
the  relationship  bewveen  the  theoretical  and  experimental  results.  The  efforts  at  Virginia  Tech  and  Florida 
Atlantic  arc  the  subject  of  a  separate  but  coordinated  project  with  ONR  HQ.  The  overall  objective  of  this 
combined  effort  is  to  evaluate  and  test  the  hypothesis  that  the  soft  owl  down  can  suppress  the  unsteady 
surface  pressure  field  on  the  scale  of  the  roughness,  either  directly  or  indirectly  through  the  modification 
of  the  turbulent  structures  of  the  boundary  layer,  and  control  not  only  the  source  of  roughness  noise  but 
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possibly  also  that  of  trailing  edge  noise  and  surface  vibration. 

At  Lehigh  University,  this  research  effort  has  defined  2  months  of  summer  work  for  Justin  Jaworski, 
which  has  resulted  in  2  conference  papers  and  1  patent  application  produced  in  collaboration  with  the 
international  research  team. 

APPROACH 

Analytical  models  for  the  suppression  and  generation  of  noise  by  the  owl  down  material  are  developed 
and  solved  to  gain  insight  into  the  essential  physical  mechanisms.  For  vortex  noise  generated  by  an 
inclined  fiber,  progress  is  made  using  conformal  mapping  techniques  coupled  with  compact  aeroacoustic 
Green’s  functions  to  find  the  far-field  sound  due  to  vortex  motion.  A  separate  analysis  concerning  noise 
suppression  by  the  presence  of  a  porous  canopy  involves  the  numerical  solution  of  the  two-dimensional 
Rayleigh  equation  using  a  standard  Newton-Raphson  numerical  routine.  Both  theoretical  analyses  identify 
the  dimensionless  groups  to  furnish  comparison  with  current  and/or  future  experiments  and  computations. 

WORK  COMPLETED 

Two  theoretical  models  are  developed  to  support  the  collaborative  investigation  of  flexible  roughness 
noise  carried  out  with  Cambridge,  Virginia  Tech,  and  Florida  Atlantic  universities.  First,  the  problem  of 
vortex  sound  produced  by  an  inclined  fiber  in  a  mean  flow  was  solved  via  conformal  mapping  techniques, 
whereby  the  computed  vortex  path  generates  a  time-dependent  signal  in  the  acoustic  far-field.  This  model 
furnishes  a  theoretical  framework  to  understand  the  noise  generation  by  owl  down  in  its  physical  limit  as 
a  sparse  collection  of  hydrodynamically-uncoupled  fibers,  enabling  further  work  to  take  into  account 
dynamic  motions  of  the  fiber  and  finite  distances  between  fibers  for  more  realistic  aeroacoustic  models  of 
its  turbulence  noise  generation. 

Experimental  measurements  made  at  Virginia  Tech  indicate  that  a  porous  canopy  can  reduce  surface 
pressure  amplitudes  by  up  to  30dB  over  a  finite  range  of  frequencies  that  depends  linearly  on  the 
background  flow  speed.  To  examine  the  hypothesis  that  the  porous  layer  creates  an  unstable  wave  near 
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the  canopy  by  inflection  of  the  boundary  layer  profile,  a  inflectional  mixing-layer  model  is  formulated 
and  solved  by  numerical  solution  of  the  Rayleigh  equation.  Numerical  solutions  indicate  that  the 
boundary-layer  inflection  point  is  pushed  upward  towards  the  canopy,  suggesting  that  a  pure  mixing-layer 
model  without  boundary-layer  (no-slip)  details  may  be  sufficient  to  describe  the  most  relevant  physics. 
The  semi-empirical  model  predictions  for  the  porous  canopy  arc  in  excellent  agreement  with  low- 
frequency  noise  suppression  trends.  However,  further  work  is  necessary  to  include  the  proper  physics  into 
the  model  to  improve  the  theoretical-experimental  agreement  at  high  frequencies. 

RESULTS 

1.  Vortex  Sound  from  Isolated  Rigid  Owl  Down  Fibers  Inclined  in  a  Mean  Flow 

The  noise  generation  due  to  the  interaction  of  the  boundary-layer  turbulence  with  the  downy  material  is 
presently  idealized  as  a  vortex  sound  problem,  where  it  is  assumed  that  the  downy  material  is  sufficiently 
sparse  such  that  its  fibers  are  hydrodynamically  uncoupled,  and  it  is  further  assumed  that  these  fibers  are 
acoustically  compact.  The  acoustical  compactness  of  the  fiber  enables  the  use  of  the  two-dimensional 
Green’s  function  derived  by  Howe  [3,  p.  138]  to  relate  the  vortex  path  near  the  fiber  to  the  resulting  sound 
radiated  to  the  acoustic  far-field. 


The  mathematical  analysis  for  sound  due  to  a  vortex  passing  round  an  inclined  fiber  follows  the  procedure 
outlined  by  Howe  [3,  pp.  193-194]  for  a  vortex  passing  over  a  rigid,  upright  spoiler,  which  includes  the 
effects  of  a  background  mean  flow  on  the  acoustic  signal.  We  extend  this  work  to  allow  for  arbitrary 
angles  of  inclination  with  (0,1),  as  illustrated  in  Fig.  1,  where  the  results  of  Howe  are  recovered  for 
d=\ll .  The  mapping  from  the  z-plane  to  the  ^  -plane  is  facilitated  by  using  an  elementary  translation  of 
the  origin  in  the  z-plane  to  the  end  of  the  spoiler,  say,  Zj  =z  +  he~'''^ ,  and  then  relating  the  -plane  to  the  ^ 
-plane  by  the  Schwarz-Christoffel  transformation  [4,5]. 
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Figure  1.  Vortex  noise  from  an  inclined  fiber.  Conformal  mapping  relates  the  physical  z-plane  to  the  f  -plane 
where  the  problem  is  solved.  The  impermeable  wall  condition  is  satisfied  automatically  by  adding  an 
appropriate  image  vortex  to  each  plane. 


The  constants  K  and  b  are  determined  by  linking  the  maps  between  the  z-  and  ^-planes.  For  example,  as 
the  same  as  for  when  ^->0.  Applying  these  conditions  to  Eq.  (1),  and  noting  that  a 
change  of  variable  to  transforms  the  integral  into  a  beta  function  B(a,/?)  [6],  these  constants  arc 

found  to  be 


b  =  \ld. 


(2) 


K  = 


hslnjzd  ^  h 


(3) 
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where  the  incomplete  beta  function  is  defined  by  [6]  , 


dr 


B{a,b) 


(4) 


Note  that  f{6)  =  \  in  Eq.  (3)  in  the  special  case  of  an  upright  fiber,  8=il2 .  The  mapping  relationship 
between  the  z-  and  f -planes  is  completed  by  equating  the  background  mean  flows  far  from  the  fiber, 
which  yields  V  KU  . 

The  acoustic  pressure  due  to  vortex  motion  requires  knowledge  of  the  vortex  path  and  the  streamwise 
component  of  the  Kirchhoff  vector  [3,  pp.  52,  192].  The  established  nonlinear  equations  of  motion  for  the 
vortex  at  ZoW  =  -^oW  +  ^>oWare  [3,  p.l07] 
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where 

|J||=-^(?(?-i)(?-i/0))-'r(zo),  (6) 

and 
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can  be  found  from  Eq.  (1)  using  the  chain  rule  or  via  direct  integration,  respectively;  the  primes  denote 

f 

(  )  ^dldz.  In  the  special  case  of  0  =  1/2 ,  Eqs.  (5-7)  can  be  combined  to  yield  equations  of  motion  that 
are  explicit  in  tenns  of  the  vortex  position,  Zq  .  However,  in  general,  Eq.  (5)  depends  implicitly  on  vortex 
position  through  f(z) ,  which  requires  the  solution  of  Eq.  (7)  for  a  given  z.  These  vortex  equations  of 
motion  (7)  can  be  expressed  in  dimensionless  form. 
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using  the  following  dimensionless  groups: 


Z  = 


(9) 


The  component  of  the  Kirchhoff  vector  in  the  streamwise  direction  is  simply  the  real  part  of  the  complex 
velocity  potential  for  unit  speed  ( t/  =  1 )  in  the  z-plane,  i.e.  Re(A'^) .  Howe  [3,  p.  193]  defines  a  function, 


such  that  the  acoustic  pressure  in  the  far-field  due  to  the  vortex  motion  is 
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where  the  second  term  in  the  integrand  implies  the  function  dependence 

A7 

— ^([7)- A“)  =  w([r]- A“)  +  /v([ri-A“)  and  f  =  fri-A^,  Equation  (11)  measures  the  acoustic  pressure  at 
dT 

far-field  location  x ,  which  is  inclined  at  an  angle  0  measured  clockwise  from  the>^-direction.  In  Eq.  (1 1), 
Pq  is  the  stationary  fluid  density,  c^is  the  isentropic  speed  of  sound,  the  Mach  number  is  defined  as 
M  =  V/cg ,  and  [7]  is  the  nondimens ional  retarded  time. 


The  vortex  path  and  acoustic  pressure  are  evaluated  numerically  using  Eqs.  (5)  and  (11),  respectively.  The 
initial  location  of  the  vortex  is  sufficiently  far  upstream  that  its  initial  motion  is  undisturbed  by  the  fiber; 
in  this  study,  the  initial  vortex  position  is  set  to  Z  =  -10  + 0.75/ ,  Figure  2  illustrates  these  vortex  paths  for 
0  =  1/3,  1/2,  and  2/3  and  in  the  cases  with  (e  =  l)  and  without  (e— 0)  a  background  mean  flow.  Curiously, 
the  vortex  path  is  pushed  upwards  by  the  fiber  at  greater  distances  than  the  respective  no-flow  cases 
regardless  of  the  fiber  angle.  Nonetheless,  the  effect  of  a  mean  flow  is  to  draw  the  vortex  path  nearer  to 
the  edge  of  the  fiber  as  it  passes  over,  resulting  in  a  louder  sound  seen  in  comparison  of  Figs.  3  and  4. 
Figures  3  and  4  show  the  reduced  acoustic  pressure  from  Eq.  (11)  as  a  function  of  the  non-dimensional 
retarded  time,  where  [7]=0  refers  to  the  instant  where  the  point  vortex  crosses  the  line  extending  from  the 
end  of  the  fiber.  The  cases  for  0  =  1/2  replicate  the  existing  results  in  the  literature  [3,  7]  and  verify  the 
new  approach  presented  here  that  is  able  to  handle  arbitrary  angles  of  the  fiber. 

A  comparison  of  the  cases  0=1/3  and  2/3  in  Fig.  4  indicates  that  leaning  the  fiber  against  the  flow 
produces  a  greater  positive  pressure  as  the  vortex  nears  the  end  of  the  fiber,  and  leaning  the  fiber  with  the 
flow  reduces  the  positive  pressure  peak  but  increases  the  magnitude  of  the  negative  peak,  as  one  would 
expect  intuitively.  However,  Fig.  3  shows  surprisingly  that  the  same  trend  exists  without  a  background 
flow,  indicating  that  the  presence  of  a  mean  flow  is  not  responsible  for  the  relative  shifts  of  the  pressure 
peaks,  although  a  mean  flow  certain  increases  the  magnitudes  of  these  peaks. 
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Figure  2.  Vortex  paths  round  a  rigid  inclined  fiber.  The  self-induced  vortex  paths  show  the  influence  of  the 
fiber  inclination  and  background  mean  flow  on  the  vortex  trajectory.  6  =1/3  (top),  1/2  (middle),  2/3  (bottom). 
The  trajectories  for  £  =  0  and  e  -  1  are  indicated  by  the  solid  and  dashed  lines,  respectively. 
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Figure  3.  Vortex  noise  signature  without  background  flow  (f  =0).  Results  for  0=  1/3,  2/3  show  that  the 
inclination  of  the  rigid  fiber  towards  the  incoming  vortex  increases  the  acoustic  pressure  values  at  the  peaks. 
The  present  methodology  is  validated  by  the  results  of  Howe  [14]  for  6  =  1/2.  The  retarded  time  is  taken  with 
respect  to  the  instant  when  the  vortex  passes  through  the  line  extending  from  the  fiber. 


Figure  4.  Vortex  noise  signature  with  background  flow  (e  =  l).  The  background  flow  greatly  enhances  both 
the  magnitude  and  approximately  doubles  the  frequency  of  the  vortex  noise  signature.  The  previous  work  by 
Howe  [1 4]  is  confirmed  for  6  -  1/2.  The  retarded  time  is  taken  with  respect  to  the  instant  when  the  vortex  passes 
through  the  line  extending  from  the  fiber. 


8 


2.  Canopy  Mixing-Layer  Instability  and  Noise  Suppression 

A  tentative  model  is  proposed  to  explain  features  of  the  reduction  in  unsteady  surface  pressure  levels 
which  result  from  introduction  of  the  fabric  covering.  First,  Fig.  5  plots  the  decibel  reduction  in  the 
experimental  surface  pressure  spectrum  between  the  datum  uncovered  case  and  the  case  with  the  covering 
in  place,  for  several  different  wall-jet  flow  speeds.  Interestingly,  the  shapes  of  the  curves  are  remarkably 
similar  to  each  other.  In  each  case  there  is  a  definite  maximum  reduction,  whose  amplitude  is  inversely 
proportional  to  I/jet  and  which  occurs  at  a  frequency  which  is  proportional  to  [/jet,  i.e.  at  constant  Strouhal 
number.  Moreover,  he  canopy  is  only  effective  in  reducing  the  surface  pressure  spectrum  over  a  finite 
range  of  frequencies  /cufon  <  f  <  /cut-off  cut-on  and  cut-off  frequencies,  /(-ut-on 

increasingly  linearly  with 

Our  inspiration  for  modeling  this  effect  comes  from  studies  of  air  flow  over  forest  canopies  (see  [8,9]  for 
a  review),  where  i:  is  well  known  that  the  drag  of  the  canopy  retards  the  lower  portion  of  the  atmospheric 
boundary  layer,  leading  to  an  inflectional  profile  which  undergoes  a  classical  mixing-layer  instability  [10] 
-  the  so-called  canopy  instability.  We  hypothesize  that  the  introduction  of  the  fabric  canopy  has  the  same 
effect  on  the  wall-jet  boundary  layer. 


The  next  step  in  our  argument  is  to  note  that  shear  instability  leads  to  exponential  growth  in  the 
streamwise  direction,  but  exponential  decay  in  the  direction  transverse  to  the  shear  layer.  This  effect  can 
be  seen  most  simply  in  the  case  of  Helmholtz  instability  of  a  vortex  sheet:  for  a  shear  of  strength  U  at 
frequency  o),  the  Kelvin-Helmholtz  wave  has  exponential  spatial  growth  rate  (o/U  and  decays  with 
transverse  distance  at  exponential  decay  rate  0)/U.  The  vortex  sheet  is  of  course  unstable  at  all 
frequencies  and  its  spatial  growth  rate  is  unbounded.  However,  a  genuine,  finite-thickness,  shear  layer 
will  typically  only  be  spatially  unstable  over  a  finite  frequency  range,  and  will  typically  have  a  definite 
frequency  at  which  a  maximum  growth  rate  (and  therefore  by  analogy  with  Kelvin-Helmholtz  a 
maximum  transverse  decay  rate)  occurs.  Moreover,  the  characteristic  frequencies  of  the  shear  layer  will 
often  occur  at  approximately  constant  Strouhal  number.  These  features  are  very  much  in  tune  with  the 
features  of  the  experimental  results  described  in  the  first  paragraph.  Our  hypothesis  for  the  effect  of  the 
canopy  is  therefore  as  follows:  the  canopy  introduces  an  inflectional  profile,  which  leads  to  instability 
waves  (as  nascent  turbulent  eddies),  which,  due  to  the  presence  of  the  canopy,  are  lifted  higher  above  the 
solid  boundary  than  they  would  have  been  without  the  canopy  present.  The  pressure  footprint  of  these 
eddies  on  the  solid  surface  is  therefore  (substantially  reduced),  by  both  the  transverse  exponential  decay 
effect  of  the  instability  wave  and  by  their  greater  stand-off  from  the  wall. 


In  order  to  investigate  our  hypothesis  further,  we  now  consider  the  stability  of  a  simple  inflectional 
boundary  layer  above  a  wall,  and  compute  the  spatial  growth  rates  of  the  system  for  varying  frequency. 
Specifically,  we  investigate  the  spatial  instability  of  the  piecewise  linear  boundary  layer,  with  a  base  flow 
which  approximates  a  general  parallel  boundary  layer  using  straight  line-segments.  Velocities  and  lengths 
are  non-dimensionalizcd  in  terms  of  the  free-stream  velocity  and  the  boundary  layer  thickness  h.  In 
terms  of  the  wall-r.ormal  coordinate  y,  the  boundary  layer  velocity  profile  is  written  as 


U(y)  =  { 


5 

{l-h)y  +  b-5 

1-5 

1 


0  <  y  <  5, 
5  <  y  <  1, 

y  >  1. 


Here  0  <  5  <  1  controls  the  wall-normal  location  at  which  the  line  segments  change  their  slope,  and 
0  <  h  <  1  determines  the  gradient  of  the  velocity  profile  in  the  two  regions  beneath  y  =  1.  An  example 
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profile  is  given  in  Fig.  6.  The  condition  for  the  profile  to  be  inflectional,  and  therefore  inviscidly  unstable, 
is  b  <  S.  This  system  is  a  generalization  of  the  case  treated  in  [1 1],  which  corresponds  to  S  =  1/2  in  our 
notation. 


Figure  5.  Experimental  surface  pressure  level  reductions  due  to  canopy  introduction. 


b 

Figure  6.  Schematic  of  the  piecewise  linear  boundary  layer. 


To  determine  the  spatial  stability  of  the  base  flow  U(y)  we  first  assume  wavelike,  two-dimensional 
velocit}^  and  pressure  perturbations  of  the  form 

[u.v,p](,x.y,t)  = 


10 


where  x  is  the  streamwise  coordinate,  o)  G  M  is  the  temporal  frequency  and  fc  6  C  is  the  streamwise 
wavenumber.  Thus  spatially  unstable  modes  have  Im{/c}  <  0.  Perturbations  of  this  form  satisfy  the 
Rayleigh  equation 


{(X)  —  kU)(D^  —  —  U"v  =  0, 

where  D  =  d/ dy.  Since  U (y)  is  linear  in  y,  and  since  /c  G  C,  the  Rayleigh  equation  simplifies  to 


(iD^-k^)v  =  0. 


This  simplification  allows  us  to  determine  an  exact  dispersion  relation  d(aj, /c,  h,  5).  Assuming  a  wall- 
normal  velocity  perturbation  profile  which  satisfies  the  Rayleigh  equation  in  each  of  the  three  mean 
profile  regions  gives 


fAe-^y  +  Be^y  0<y  <5, 
viy)  =\ce-^y  A- De^y  5  <  y  <  1, 

I  Ee~^y  y  >  1; 

for  Re{/c}  >  0.  We  can  solve  for  d((ji),k,h,S')  by  elimination  of  the  constants  A,B,C,D,E.  These 
constants  are  eliminated  by  imposing  zero  wall-normal  flow  at  the  lower  wall,  and  jump  conditions  at  5 
and  1  which  impose  continuity  of  pressure  and  wall-normal  velocity.  Using  IJ  to  denote  jump 
conditions,  we  have 


i;(0)  =  0, 

Kcl)  —  kU)Dv  +  kU'v}  =  0, 


In  this  manner,  we  derive  the  four  parameter,  nonlinear  dispersion  relation 

d!(cu,  k,  b,  S')  =  acj^  +  Po)  +  y, 
where,  using  ^  =  coth(/cd),  we  have  defined 


^  =  (1  +  0 


a  =  2(l  +  a 


1-h 
1  -  d 


-  2k{l  +  h) 


2(b  -  S) 
^5(1-5) 


r  =  (6«1  +  f )  +  (zk  -  +  (^)  -  M(1  -  o) 


1-b 


1-b^ 


S-b 


5(1  -  5) 


5(1  -  5) 


Solutions  k  of  5(0),  k,  b,5)  =  0  are  found  via  Newton  iteration,  where  the  derivative  dd/dk  is  computed 
at  each  step  from  an  analytic  expression. 


The  theoretical  framework  presented  above  is  now  used  to  attempt  to  reproduce,  as  far  as  possible,  the 
experimental  results.  In  order  to  do  this,  we  will  make  sensible  choices  for  some  of  the  parameters  in  our 
model  which  can  be  estimated  (such  as  the  jet  speed  and  the  canopy  height).  However,  other  parameters 
(such  as  the  effective  strength  and  position  of  the  inflection  point)  are  harder  to  estimate,  and  instead  we 
will  try  to  tune  the  free  parameters  b  and  5  in  order  to  try  and  fit  the  results.  We  are  therefore  not 
attempting,  at  this  stage,  to  make  a  complete  quantitative  prediction  without  the  need  for  any  fitting,  but 
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rather  we  wish  to  determine  whether  or  not  the  canopy  instability  model  is  capable  of  recovering  the 
features  seen  in  experiment. 

To  facilitate  experimental  comparison,  the  boundary  layer  thickness  in  our  model  is  taken  to  be  the 
canopy  height  h  =  =  4mm.  The  velocity  above  the  canopy  is  observed  experimentally  to  be  about  one 

third  of  the  total  jet  velocity,  so  accordingly  our  free-stream  velocity  is  chosen  as  =  t/jet/3. 
Frequencies  /  in  Hertz  can  then  be  mapped  to  non-dimensional  frequencies  using  the  relation 

6nfhc 

0)  = - . 

t/jet 

For  each  the  parameters  b  and  S  are  tuned  by  matching  the  cut-on,  cut-on'  maximum,  Wmaxi 
frequencies  of  the  instability  growth  rate 

Im{/c(aJeut-on)}  =  0,  lm{k(^0  <  Wcut-on)}  >  0, 

minlm{/c((ij)}  =  Im  {/cCWfuax)}- 

with  the  experimental  values  and  w^ax-  optimize  over  b  and  S  to  satisfy  the  relative  error 

constraint 


Results  of  the  optimization  arc  shown  in  Fig.  7,  where  it  can  be  seen  that  the  optimal  b  lie  in  the  range 
0.825  <  b  <  0.855  and  the  optimal  S  are  scattered  between  0.905  <  S  <  0.945.  We  plot  the  growth 
rates  in  Fig.  8. 
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The  pressure  eigenfunctions  for  each  unstable  eigenvelocity  can  be  computed  using  the  relation 

Examples  for  cases  close  to  the  cut-on  and  maximum  frequencies  are  plotted  in  Fig.  9.  Note  in  Fig.  9(a) 
that  the  pressure-amplitude  of  the  instability  wave  changes  little  between  just  above  the  canopy  and  the 
wall,  in  the  cut-on  case.  However,  there  is  significant  attenuation  in  the  least  stable  case,  seen  in  Fig.  9(b), 
in  support  of  our  assertion  that  the  presence  of  the  canopy  shields  the  wall  from  the  presence  the  pressure 
fluctuations  in  the  upper  boundary  layer. 

In  order  to  compare  more  directly  with  the  SPL  attenuation  data,  we  suppose  that  the  change  in  SPL  is 
related  to  the  linear  instability  in  the  following  manner 

ASPL  =  201ogio(e“‘"’^''^‘^''), 

where  we  define  dj^  to  be  the  pressure  instability  length  scale,  i.e.  the  height  above  the  wall  where  the 
least  stable  linear  wave  is  located  in  order  to  fix  the  level  of  attenuation  observed  in  the  experiment.  For 
each  jet  speed,  we  determine  using  the  maximum  instability  growth  rate  and  the  maximum  SPL 
fluctuation 


dh  — 


-log  10 


flO  20  j. 


In  Fig.  6(a)  we  plot  the  value  d^  takes  at  each  jet  speed,  and  in  Fig.  6(b)  we  contrast  our  ASPL  estimate 
with  the  experimental  quantities. 
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Figure  9.  Pressure  eigenfunctions  for  =  40  (a)  o>  =  1, 2  (b)  to  =  8. 05. 
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Figure  10.  (a)  values  of  the  pressure  instability  length  and  (b)  our  ASPL  estimate  (solid)  against  the 

experimentally  observed  values  (dashed). 


What  becomes  apparent  in  Fig.  10(b)  is  that  the  canopy-instability  model  is  capable,  for  suitably  chosen 
parameter  values,  of  capturing  aspects  of  the  canopy  attenuation,  specifically  the  lower-frequency  cut-on 
behavior  and  the  maximum  attenuation.  However,  although  the  canopy-instability  model  does  indeed 
predict  a  maximum  frequency  (and  captures  the  right  trend  with  varying  speed),  the  predicted  value  of  the 
maximum  frequency  differs  significantly  from  the  experimental  value.  Indeed,  it  is  clear  that  the 
experimental  attenuation  curves  differ  in  shape  from  the  canopy-instability  curves  above  a  certain  speed- 
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dependent  frequency  (note  the  ‘shoulder’  in  each  experimental  curve),  which  suggests  that  some  other 
mechanism  takes  over  at  higher  frequency.  Elucidation  of  this  other  mechanism  is  a  matter  for  further 
research. 

IMPACT/APPLICATIONS 

The  theoretical  models  developed  through  this  project  offer  insight  into  the  noise  suppression  behavior  of 
individual  fibers  and  porous  canopies,  which  may  impact  further  research  into  turbulence  noise  generation 
and  its  interactions  with  more  sophisticated  poroelastic  wall  structures. 

The  international  research  collaboration  led  to  the  development  of  streamwise-oriented  surface  structures 
termed  “finlets,”  which  have  not  been  included  in  this  report  due  to  its  patent  application  status.  Our  flnlet 
designs  have  been  demonstrated  experimentally  to  yield  trailing-edge  noise  suppression  levels  of  up  to 
lOdB  over  a  broad  range  of  angles  of  attack  with  minimal  impact  on  the  aerodynamic  lift  and  drag 
characteristics.  The  physics  mechanisms  underpinning  the  flnlet  design  are  not  completely  understood 
and  remain  the  subject  of  continued  theoretical  and  experimental  investigations. 
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